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We investigate the effect of early chemical freeze-out on hydrodynamic flow and HBT 
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1. INTRODUCTION 

From the particle ratios and the transverse momentum spectra of hadrons in relativis- 
tic heavy ion collisions, one can obtain chemical freeze-out temperature T ch and thermal 
freeze-out temperature T th through the analyses of statistical models These temper- 
atures are usually different from each other at the AGS, SPS and RHIC energies, i.e., 
T ch ~ 160-200 MeV while T th ~ 100-140 MeV. This provides the following picture of 
space-time evolution of hadronic matter: The equilibrated system produced in the early 
stage of collisions first undergoes the chemical freeze-out where the observed particle ratios 
are fixed and then goes through the thermal freeze-out where the momentum distribution 
is fixed. However, one assumes T ch = T th in the conventional hydrodynamic models. As a 
result, one cannot describe the particle ratios and spectra simultaneously. As a first trial, 
we incorporate the early chemical freeze-out into the hadronic equation of state (EOS) 
for /i B = case || and investigate its effects on hydrodynamic evolution such as radial 
flow, elliptic flow and HBT radii by using a fully three-dimensional hydrodynamic model 
(3D hydro) §. 

2. MODEL EOS 

We construct two models EOS to compare the space-time evolution of fluids. Both 
models describe the first order phase transition between the QGP phase and the hadron 
phase at T c = 170 MeV. Here we take the baryonic chemical potential ^ = for sim- 
plicity. The QGP phase is composed of massless free u, d, s quarks and gluons and is 
common to both models. We employ the following two models for the hadron phase. One 
is the conventional resonance gas model in which chemical equilibrium is always assumed 
(the model CE). We include all baryons and mesons up to the mass of A(1232). The 
other model describes the partial chemical equilibrium (the model PCE) in the hadron 
phase m. Whether an interaction is chemically equilibrated in the hadron phase depends 
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Figure 1. Transverse momentum distributions of negative charged hadrons in Au + 
Au central collisions at ^s^n — 130 GeV for the model CE (left) and for the model 
PCE (right). Solid, dotted and dashed lines correspond to T th =140, 120 and 100 MeV, 
respectively. 



on its time scale. Elastic scattering through a resonance particle, such as — > p — ► ttty, 
ttK —>■ K* —>■ tcK and ttN —>■ A — ► ttN, has large a cross section in the hadronic 
medium and frequently happens within the lifetime of fluids. Hence these interactions 
are equilibrated in the model PCE. On the other hand, particles which decay through 
the weak interaction have much longer lifetimes than the typical time scale of the fluid 
(~ 10 fm/c). So those particles are assumed to be chemically frozen in the model PCE. 
Below T ch , we can introduce chemical potential /Zj associated with the conserved parti- 
cle numbers. The ratio of the conserved number including contribution from resonance 
decays iVj(= iVj + J2j^idj^ixNj) to the entropy S has a constant value along an adi- 
abatic path. Thus we obtain the chemical potential as a function of temperature from 
the equation rij(T, ^)/s(T, = nj(T ch ,/ij = 0)/s(T ch , {/^} = 0), where n and s are, 
respectively, the particle number density and the entropy density. We regard it, K, rj, N, 
A, S and its anti-particles as "stable" particles and chemical potentials for all hadrons are 
represented by chemical potentials associated with these stable hadrons, e.g., \i p = 2/x T , 
\ik* = P-k + (J>k, A* A = /Wtt + and so on. We note that the chemical potential of 
anti-nucleons is equal to that of nucleons in this model: = A*at- 

3. RESULTS 

We choose initial parameters so as to reproduce the pseudorapidity distribution of 
charged hadrons for central (0-6%) and non-central (35-45%) Au + Au collisions at 
^snn = 130 GeV |J. The resultant maximum value of the initial energy density in 
collisions with vanishing impact parameter is 35 GeV/fm 3 at the initial time To = 0.6 
fm/c. For collisions with finite impact parameters, the initial energy density is assumed 
to scale with the number of binary collisions. For details, see Ref. ||. The transverse 
momentum spectra of negative charged hadrons observed by the STAR Collaboration |7j 
are compared with our results in Fig. 0. It is found that the pt slope is insensitive to 
the thermal freeze-out temperature T th for the model PCE while the p t slope from the 
model CE depends on T th . This is a consequence of the suppression of radial flow due to 
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Figure 2. 1*2(77) of charged hadrons in Au + Au central collisions at v /J/v r /v = 130 GeV for 
the model CE (dashed) and for the model PCE (solid). 



the early chemical freeze-out 0. We reproduce the p t spectra below 1.0-1.5 GeV/c and 
underestimate in higher p t region in the model PCE analysis. This indicates that the hard 
contribution becomes important above 1.5-2.0 GeV/c. It should be noted that contribu- 
tion from non-thermalized high p t partons with an appropriate amount of energy loss can 
fill in the difference between hydro results and the experimental data in the moderate high 
p t region (2 < p t < 5 GeV/c) [jlOl . We next show in Fig. |]the pseudorapidity dependence 



of the second Fourier coefficient of azimuthal distribution V2(rj). Here we choose T th = 140 
MeV for both models. Experimental data by PHOBOS || and STAR || are compared 
with our results. We see, as well as the radial flow, the elliptic flow is suppressed by 
taking into account the early chemical freeze-out. The shape of 1*2(77) is sensitive to the 
initial parameterization of energy density in the longitudinal direction. Bumps in forward 
and backward rapidity regions result from the crescent shape of initial energy density in 
transverse plane at space-time rapidity | T] s |~ 3-4 j4j. We might check whether other 
initial parameters can lead us to reproduce dN/drj and 7^(77) simultaneously. Finally, let 
us see the effect of the early chemical freeze-out on the HBT radii. Figure |3] shows -R S idc, 
.Rout; -Riong and -Rout /-Rsido as functions of the transverse momentum of negative pion pair 
Kt- We compare our results with experimental data by STAR JTTJ] and PHENIX [12 



-R out and -Ri on g are reduced by the early chemical freeze-out. This mainly comes from the 
reduction of the lifetime of a fluid ||. We reproduce the K T dependence of -Ri ong by the 
model PCE, so the longitudinal dynamics can be correctly described by a hydrodynamic 
model with the finite longitudinal size and the chemical non-equilibrium property. On 
the other hand, we cannot reproduce the transverse radii even in this model. Therefore 
the early chemical freeze-out is not the complete interpretation of the RHIC HBT puzzle. 

4. SUMMARY 

We investigated the effect of early chemical freeze-out on collective flow and HBT radii 
at the RHIC energy with a full 3D hydrodynamic model. It is found that radial flow, 
elliptic flow, .Rout and -Ri ong are reduced by considering the early chemical freeze-out. In 
this analysis, we took /ib = as a first trial calculation. It might be interesting to con- 
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Figure 3. HBT radii in Au + Au central collisions at ^/snn = 130 GeV. The solid 
(dashed) lines correspond to the results from the model PCE (CE). 



struct a model PCE with /xb ^ and to discuss particle ratios and spectra simultaneously 
at various collision energies. This will be discussed elsewhere. 
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